Multi-critical behavior of Z 2 x 0(2) Gross-Neveu- Yukawa theory in graphene 
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Multi-critical behavior of interacting fermions in graphene's honeycomb lattice is presented. In particular, 
we considered the spin triplet insulating orders, where the spin rotational symmetry of the order parameter is 
explicitly broken. By casting the problem in terms of Gross-Neveu- Yukawa theory we show that such symmetry- 
breaking terms are irrelevant near the metal-insulator critical point. A finite Yukawa coupling among bosons 
and fermions improves the stability of such critical point against the symmetry-breaking perturbations. Physical 
sources of such symmetiy-breaking terms are pointed out. Critical exponents are calculated near the transitions 
as well. 
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The pseudo relativistic nature of low energy quasi parti- 
cles opened a new frontier in condensed matter physics with 
the discovery of graphene, a monolayer of graphite. 1 Even 
though in its pristine state graphene appears to be a semi- 
metal, fermions may also find themselves in various insulating 
phases as well. Depending on the relative strength of the dif- 
ferent finite range components of the Coulomb interaction, the 
ground state can lack a plethora of discreet or continuous sym- 
metries. For instance, a sufficiently strong onsite Hubbard in- 
teraction (U) orients the spin at each site in an opposite direc- 
tion from its neighboring ones, whereas the average electron 
density acquires a staggered pattern at large nearest-neighbor 
Coulomb repulsion (Vi)£^ Yet another insulating phase may 
result from a strong next nearest neighbor Coulomb repulsion, 
which can induce a gapped insulating phase with finite cir- 
culating currents between the sites on the same sublattice. 5,6 
This state is named the quantum anomalous Hall (QAH) state 
and breaks time reversal symmetry (TRS) upon acquiring a 
finite expectation value in the ordered state. However, fluc- 
tuations preempt appearance of a QAH state and stabilize the 
quamtum spin Hall (QSH) insulator. The QSH state breaks 
the TRS only for each spin component. 7,8 These transitions 
out of the symmetric semi metallic phase into the gapped in- 
sulating phases are believed to be continuous and belong to 
the Gross-Neveu universality class. 3 ' 4,9 Recently it has been 
argued that Kekule bond density wave order may also appear 
when the nearest neighbor and second nearest neighbor inter- 
actions are comparable. 10 If the net interaction has an attrac- 
tive component, fermions may find various superconducting 
states to condense into. For example, an onsite attraction fa- 
vors a spin singlet s-wave superconducting ground stated. A 
second nearest neighbor attractive interaction leads to a spin 
triplet f-wave superconductor. 12 A spatially inhomogeneous 
superconducting ground state, which breaks the translational 
symmetry of the honeycomb lattice into Kekule pattern of 
bond order parameter is proposed as a variational ground state 
for strong attractive interactions between the electrons living 
on nearest neighbor sites. £11 

In recent studies^ 9 ., it has been argued that near the 
semimetal-insulator quantum critical points the pseudo rel- 
ativistic invariance of the non-interacting theory is restored. 
Consequently, the dynamical critical exponent (z) is exactly 
equal to unity leading to a non-critical behavior of Fermi ve- 
locity. There the nature of the quantum phase transitions in 



both spin singlet and triplet channels is considered, assum- 
ing that the interacting Hamiltonian describing the transitions 
does not break the spin rotational symmetry. However, the 
symmetry can be broken by, for example, the crystal strain or 
finite spin-orbit coupling. 7 Once we introduce an anisotropy 
in the order parameter along one spin direction, theory lose s 
the 0(3) symmetry associated with the spin rotation, and en- 
joys a reduced x 0(2) symmetry. In the present discus- 
sion, I am concerned with the relevance of the broken spin 
rotational symmetry of the order parameters near the quan- 
tum criticality. Within the framework of the one loop e— 
expansion, I found that such anisotropy in the spin degrees 
of freedom is irrelevant in the vicinity of the metal-insulator 
quantum critical point. However in a similar situation in an 
interacting bosonic system, it has been shown that the order- 
disorder transition is driven by a biconal fixed point, enjoying 
the Z2 x 0(2) symmetry. 15 In this discussion, I will show 
that the critical point associated with the disorder-order tran- 
sition exhibits full rotational symmetry. However, this may be 
a consequence of the one loop calculations. An interesting ob- 
servation is that in presence of finite Yukawa coupling among 
the Dirac fermions and the self interacting Higgs fields, such 
anisotropy acquires an additional degree of irrelevance near 
the transitions. Here the Higgs bosons are composite fields. 

Let us first consider a collection of free fermions on 
graphene's honeycomb lattice. The tight binding model for 
spin- 1/2 fermions, on graphene honeycomb lattice with only 
nearest-neighbor hopping, is defined as 
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where u a , v a are fermioninc operators on two triangular sub- 
lattices of the honeycmb lattice, and b{S, with i = 1,2,3 con- 
nect each sites on A-sublattice with its three nearest-neighbor 
ones. Keeping the Fourier modes near two inequivalent Dirac 
points, located at the corners of the Brillouin zone 2 , at K = 
(1, l/v / 3)27r/av / 3 and K' = —K, let us construct an 8— 
component Dirac fermion ^ = V &-) T , with 

Vj(q) = [u a (K + q),v a (K + q),u a (-K + q),v a (-K + q) 

(2) 

with a = ±, which corresponds to electrons spin projection 
along the z-axis. The tight-binding Hamiltonian in the low 
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energy approximation then takes the form 

H t = J2^(Q)M(q) + 0(q% (3) 

Q 

with Ho being the Dirac Hamiltonian in two dimensions 
which in our representation acquires a simple form 

Hd = <?o <8 ijo^iqi- (4) 

Here the four-component anticommuting gamma matrices be- 
long to the representation 70 = tq ® T3, 71 = t 3 ® t 2 , 72 = 
To ® rx, where To is the two componenet identity matrix and r 
are the standard Pauli matrices. 4 We also define the remaining 
two gamma matrices 73 = t\ (81x2,75 = t<z%ti- The two com- 
ponent Pauli matrices (<tq, <j) act on the spin indices. Here, 
for convenience we set the Fermi velocity vp = tav3/2 to 
unity. Before we consider the effect of electron-electron inter- 
action on the quasi-particle dispersion, it is useful to register 
the symmetries of the free Hamiltonian. The Dirac Hamilto- 
nian commutes with P = 00 (8*7375, which in our representa- 
tion stands for the generator of translation. It also commutes 
with Ik = (To ® 17175 and I uv = ctq (8 72, upon inverting 
the momentum axes q\ — > —q\ and qi — > —q^, respectively. 
These two operators, respectively correspond to the exchange 
of two non-equivalent Dirac points and the sub-lattices. 

We now consider the effect of electron-electron interaction 
on the gapless excitations spectrum in the vicinity of the Dirac 
points. Due to vanishing density of states at Fermi energy, 
the gapless excitations are robust against any small electron- 
electron interaction. However, if the interactions are suffi- 
ciently strong, the system may suffer semi-metal insulator 
transitions. In 2 + 1 dimensions there are two ways of gener- 
ating a dynamical mass at sufficiently strong interactions. The 
first order-parameter (OP) is 

4> = (4> s ,$t) = ((* f cr (870*), (* t ff(8 7o*))- (5) 

4> s (4>t) preserves (breaks) the TRS and breaks chiral C/ c (4) 
symmetry (CS) generated by {ct ,ct} (8 {h, 73, 75, 735}, 
where 735 = 17375. The TRS is defined as ^ a — > It^a, 
where I t is an anti-unitary operator defined as I t = UK, K 
being the complex conjugate. Here U is the unitary part and 
in 'graphene representation' U — 17175. In our representa- 
tion <j) s corresponds to a finite chemical potential differing in 
its sign from its neighbors, whereas <fi t to a finite Neel order- 
ing. The second order parameter 

X = (Xs,Xt) = ((* t cr (8i7i72*),(* t ff® «7i72*>), (6) 

on the other hand breaks the TRS for each spin component but 
preserves the CS. Such a correlated ground state supports an 
intra sub-lattice current, introduced by Haldane 5 , circulating 
in opposite directions on two sub-lattices. 

It is worth mentioning that the singlet components of the 
OPs (xs, 4>s) break the Ising like symmetry among the sub- 
lattices, whereas the triplet OPs (xt,<fit), additionally break 
the full spin rotational symmetry. Therefore the quantum crit- 
ical points (QCPs) corresponding to the development of \s 



and (fi t belong to the same universality class. Analogously the 
QCPs associated with the generation of \t and (fit also en- 
joy the same universality class. 9 Moreover all the OPs listed 
above transform as scalars under the pseudo Lorentz transfor- 
mation. Therefore near the quantum critical points, driving 
the system out of the semi metallic phase to the ordered insu- 
lating phases the dynamical critical exponent(z) is expected to 
be unity. Nevertheless, a Loretz symmetry breaking perturba- 
tion is found to be irrelevant near the QCPs4£ 

Next we concentrate on the spin triplet insulating orders, 
namely (xt, (fit) an d at the same time set \s — <fis — 0. Let us 
start with the onsite Hubbard interaction (U) only, for simplic- 
ity. However, we introduce a non-trivial anisotropy along one 
particular spin projection. On methodological level, it appears 
interesting to study the relevance of such anisotropy near the 
semimetal-insulator quantum critical point. The quantum me- 
chanical action in the presence of anisotropic interaction along 

one particular spin projection reads as S — J^ T drdxLi nt , 
where 

Unt = *o- ®7M 5 M*+9||(* t<7 ll« ) 7o*) 2 +.9±(* t cr3®7o*) 2 , 

(7) 

where an = (<ti,(T2) and fj, = 0,1,2. fj, = stands for the 
imaginary time component. The Einstein summation conven- 
tion is assumed, but only over the repeated space time indices. 
Within the framework of the Hubbard model with only onsite 
repulsion one discovers g\\ = g± = U 

By performing the Hubbard-Stratonovich transformation 
one can rewrite the effective action corresponding to Li nt in 
d— dimension as 

S = J d d x{~[^a + g\\(fi\\{& &\ (870*) 

+ g±(fi±{^a 3 ® 70*)] + ^[(d„(fi\i) 2 + (d^) 2 } + m^f 

1=1,2 i=l,2 

(8) 

Here the fermion bilinears are coupled to the Higgs OPs, ^s. 
In the Yukawa form the theory is renormalizable in 3 + 1 di- 
mension, where both Yukawa (g^g^) and Higgs self energy 

couplings (Am , Aj_, A) are essentially dimensionless. In partic- 
ular, when g 2 = g\, concomitantly A|| = A^ = A, (fi corre- 
sponds to Ne'el order parameter. 

Next, we consider the renormalization group study of the 
coupling constants in the Yukawa theory, Eq. (H)- Here we re- 
strict ourselves only to the one loop expansion. In particular, 
we are interested in studying its flow equations in d = 4 — e 
dimensions. The coupling constants (X±,Xn,X,g±,gn) are 
dimensionless in d = 4, and hence a controlled perturba- 
tive expansion can be performed in terms of a small param- 
eter e, define above 1 —. For similar studies in a system of 
interacting bosons, readers may consult Ref. 18 . Subscribing 
to the fact that the pseudo relativistic invariance is respected 
near the QCPs 4,9 , upon integrating out the fast bosonic and 
fermionic modes within the four-momentum shell A/6 < 
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FIG. 1: Flow diagram in the Yukawa coupling space. Fully stable 



fixed point is located at gt = g ± 



2JV+1 5 



for N = 2. 



equal (= e/5 for N — 2), we found the fixed point, in the 
(An , Ax, A) plane stable from all direction, resides at 



A|| = \±_ = A = 



48 
55* 



(15) 



Therefore only the fixed point with all three Higgs and both 
the Yukawa couplings having equal strength has one unstable 
direction and thus is critical^ The unstable direction corre- 
sponds to flow of the mass of the Higgs field, m\\ = m± = m 
and its flow reads as 



dm 2 



= 2m z 



-Am 
6 



2Ng 2 m 2 , 



(16) 



{lu 2 + k 2 ) 1 / 2 < A (b > 1), one can write down the differ- 
ential flow equations of the coupling constants as 



d\nb 



dgj 

d\nb 



egf l -(2N + 2)gf l +g 2 g 2 ± , (9) 



eg 2 ± -(2N + 3)gi + 2g 2 g 2 ± , (10) 



dX \\ 
d\nb 



eAy - (- Af + - A 2 ) - 4iVA||.gf + 2ANg\, (11) 



|^ = eX x (§Ai + \x 2 ) - ANX^gl + 24Ngi, (12) 



dX 
d\nb 



eX- (AAii + — An A_l + :;A 2 ) 



(i 



MNgfgi, 



2NX ± (gi+g 2 ) 
(13) 



where^ = S d A e g x , X x = S d A e X x , x =||,_L, and A = 
SdA e X. Sd is the surface of a d— dimensional unit sphere and 
N corresponds to the number of 4— component spinor. Thus 
for graphene N is equal to 2. Here the ultraviolet cut-off, 
A « 1/a, corresponds to the interval of energy over which the 
linear approximation of the density of states holds. 19 First we 
study the flow of the Yukawa couplings in the (e/j| , e?2_ ) pl ane - 
Since we consider the flow of the coupling constants in the 
critical hyperplane, we wish to find the fixed point which is 
stable from all directions. From the first two flow equations, 
Eq. 191101 . one finds that the fixed point in this plane stable 
from all direction is located at 
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9l 



2N + 1 



(14) 



Therefore near the quantum critical point, two Yukawa cou- 
plings, (<7 2 ,<Ti_) enjoy equal strength, Fig.Q] This feature is 
independent of the value of N. Hence for the rest of our dis- 
cussion we set g 2 = g\. Next we consider the Higgs sec- 
tor of the theory. Upon setting both the Yukawa couplings 



where g 2 — g\ — g and A|| = A^ = A = A. The mass of the 
Higgs field (m) is proportional to the temperature. Therefore, 
near the metal-insulator quantum critical point anisotropic 
coupling turned out to be irrelevant and consequently the spin 
rotational symmetry is restored. The emergence of 5*0(3) 
symmetry near the QCP happens to take place for arbitrary 
N. A renormalization group study on a closely related, but 
purely bosonic $ 4 theory (py = g± = 0), is performed over 
a most general 0(ni) © 0(^2) symmetric Landau-Ginzburg- 
Wilson Hamiltonian, involving two fields <fii and (f>2 enriched 
by 7ii and components, respectively. In an extensive five 
loop e— expansion it was found that for ri\ = 1 and 712 = 2, 
the disorder-order transition is governed by a biconal fixed 
point possessing a^x 0{2) symmetry. However, the critical 
exponents near the biconal fixed point are extremely close to 
the ones for the 5*0(3) symmetric Heisenberg fixed point. 15 
It is worth mentioning that the critical point associated with 
the order-disorder transition enjoys the 0(2) symmetry even 
though the interactions are invariant under a Z2 x Z2 sym- 
metries. On the other hand, when ri\ + ri2 > 3, the system 
finds itself in an ordered phase via decoupled fixed points. 
It is admitted that the emergence of the spin rotational sym- 
metry near the quantum critical point may be an artifact of 
the one loop calculations. However, it is worth mentioning 
that the negative eigenvalues of the stability matrix near the 
spin symmetric fixed point announced above, acquires addi- 
tional contributions of the same sign for finite Yukawa cou- 
plings. Therefore, the flow of the irrelevant trajectories to- 
ward the critical point is faster when the Higgs fields are cou- 
pled to the fermions. By taking g^ = g± = the eigen- 
values of the stability matrix at A|| = A^ = A = ||e are 
(—1, —7/11, —2/11) e, whereas those with 9|| = g± = e/5 
are (-19/5, -177/55, -137/55) e, with N — 2. Therefore 
it may be worth studying the multi-critical behavior of the 
present model beyond one loop level. 

Finally let us consider the other triplet order \t = (Wcf® 
*7i72^), that breaks the TRS for each spin component. The 
insulating state is referred as quantum spin Hall insulator 
(QSHI). Upon incorporating the fluctuations in the ordered 
phase around the saddle point, it also preempts appearance of 
quantum anomalous Hall insulator. A sufficiently strong sec- 
ond nearest neighbor repulsion can take the system into such 
an ordered phase. 5,6 The spin rotational symmetry of this or- 
der parameter is broken in the presence of a finite spin-orbit 
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coupling, which is proportional to the third component of the 
order parameter. 7 In a graphene system the spin-orbit coupling 
~ 0.01 — 0.2K, is extremely small in comparison to the finite 
range Coulomb repulsions. 21 Even though we cannot study 
the critical properties of the semimetal-insulator instability to- 
ward the QSHI phase, at d = 2 + 1 it belongs to the Gross- 
Neveu universality class of the anti-ferromagnetic instability. 
This is because the term ^[cr ® 27375] 'I' can be transformed 
into ^[ct®/^^, while keeping the free Lagrangian invariant. 4 
Therefore, in the presence of a finite spin-orbit coupling the 
order parameter Xt li yes on the surface of the S2, sphere in 
three dimensions. 

To summarize, we here considered the spin triplet insulat- 
ing order of Dirac fermions, when the spin rotational symme- 
try of the order parameter is explicitly broken. Performing the 
Hubbard-Stratonovic transformation, we present the theory as 
a Z2 x 0(2) symmetric Gross-Neveu- Yukawa theory. Within 
the framework of one loop e— expansion, where e = 4 — d, 
we found a quantum critical point that restores the symme- 
try under spin rotation and drives the system from a symmet- 
ric semimetallic phase to a ordered insulating phase. 9 Even 
though the conclusion is based on a simple one loop calcu- 
lation, we found that such QCP acquires extra stability when 
the bosonic Higgs fields are coupled to fermionic fields via 
Yukawa interactions. This result dictates that even when there 
exists anisotropy in an interacting model at the lattice scale, 
along different spin directions, it smears out near the transi- 
tion and restores the symmetry under spin rotation. Recent 
quantum Hall experiments 22 , revealed a Kosterlitz-Thouless 
scaling of the resistivity near the metal insulator transition, 
when the system is tuned to filling one-half. This observa- 
tion initiated a search for the origin of the vortex like exci- 
tations of the Dirac fermions. One candidate for the possi- 
ble order parameter with the requisite ?7(1) structure is the 
Kekule bond-density-wave, which is favored by the electron- 
phonon coupling. 23,24 Another possibility is the antiferromag- 
netic ordeF 2 ^ 2 ^, favored by onsite Hubbard interaction, pro- 



jected onto an easy plane by Zeeman coupling. In the pres- 
ence of a magnetic field, the relativistically invariant band, 
collapses onto set of Landau levels (LL) and existence of a 
half-filled LL drives the system through metal-insulator tran- 
sition even at infinitesimal interactions. At zero field critical- 
ly gap shows a perfect square root dependence on the mag- 
netic field. If the Hubbard interaction is not too far from its 
critical (zero field) strength for insulation 19,27 , scaling (sublin- 
ear) of the interaction induced gap is essentially determined 
by the zero field critical properties. The magnetic field then 
plays the role of finite size length scale. Such scaling function 
has been computed field theoretically, as well as numerically 
in the Hartree limit 25,28 . In this work, I showed that the Neel 
order is unlikely to be restricted in the easy plane, within the 
frame work of the anisotropic Hubbard model. It is admitted 
that the source of anisotropy in the Hubbard model is currently 
unknown. However, a finite Zeeman coupling immediately re- 
stricts the Neel order in the easy plane and may support vortex 
excitations:^^ 

Finally one can compute the critical exponents near the crit- 
ical point which restores the spin rotational symmetry. These 
have been also previously computed in Ref. 9 . Namely, the 
correlation length exponent can be calculated from Eq. |[T6l , 
yielding v = 5 + Moreover one can compute the bosonic 
anomalous dimension, which is T]b = f e, and the fermionic 
anomalous dimension = as well. The fermionic 
anomalous dimension determines the behavior of the fermion 
propagator G7 ~ (a; 2 + fc 2 )( 1 ^^)/ 2 , near the critical point 
as one approaches the QCP from the semimetallic side of the 
transition. 
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